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We investigate the reconstruction of a Fermi surface, which is called a Lifshitz transition, in magnetically ordered 
phases of the periodic Anderson model on a square lattice with a finite Coulomb interaction between / electrons. We 
apply the variational Monte Carlo method to the model by using the Gutzwiller wavefunctions for the paramagnetic, 
antiferromagnetic, ferromagnetic, and charge-density-wave states. We find that an antiferromagnetic phase is realized 
around half-filling and a ferromagnetic phase is realized when the system is far away from half-filling. In both magnetic 
phases, Lifshitz transitions take place. By analyzing the electronic states, we conclude that the Lifshitz transitions to 
large ordered-moment states can be regarded as itinerant-localized transitions of the / electrons. 


1. Introduction 

The Fermi surface is an important ingredient for character¬ 
izing a metallic state. In general, the Fermi surface is affected 
by a phase transition such as a magnetic transition. On the 
other hand, the possibility of a phase transition described by 
a change in the Fermi surface topology itself has been pro¬ 
posed by Lifshitz.^ In recent years, such Lifshitz transitions 
have been discussed as a possible origin of some anomalies 
in heavy-fermion systems, for example, the phase transition 
between ferromagnetic phases of UGe 2 under pressure, 
YbRh 2 Si 2 under a magnetic field,and the transition be¬ 
tween the antiferromagnetic phases of CeRhi_j:COj:In 5 .^^ Re¬ 
cently, Fermi surface reconstruction in the antiferromagnetic 
phase of CeRhIns under a magnetic field has also been re¬ 
ported.'^ 

Such a possibility of the existence of a Lifshitz transition 
under a magnetic field and/or in a magnetically ordered state 
in /-electron systems has been investigated theoretically for 
a long time. Fermi surface reconstruction in an antiferromag¬ 
netic phase is found in the Kondo lattice model*^“^^ and in 
the periodic Anderson model.Under a magnetic field or in 
a ferromagnetic phase, there is a possibility of realizing a half- 
metallic state, where only one spin band has a Fermi surface. 
In the other phases, both spin bands have Fermi surfaces, and 
thus a transition to the half-metallic state from any of the other 
states inevitably accompanies a change in the Fermi surface 
topology. Indeed, such transitions to the half-metallic state 
have been found in the Kondo lattice modeP'*"^' and in the 
periodic Anderson model. 

In these theoretical studies, while the models are similar, 
the antiferromagnetic and ferromagnetic cases are treated sep¬ 
arately except for a Kondo lattice model with the explicit in¬ 
clusion of antiferromagnetic and ferromagnetic Heisenberg 
interactions.^^ For a further understanding of the Lifshitz tran¬ 
sitions, it is desirable to obtain a unified picture for both the 
magnetic cases. In addition, in the above studies on the pe¬ 
riodic Anderson model, the Coulomb interaction U between 


/ electrons is taken as t/ — > oo except in the studies of the 
transition to the half-metallic state by the slave-boson mean- 
field approximation^^’ and by a type of Gutzwiller approx¬ 
imation.^^ We also note that the Kondo lattice model is an 
effective model of the periodic Anderson model in the limit 
of t/ — > oo. Thus, it is unclear whether the Lifshitz transition 
exists or not even for a finite U beyond these approximations. 

In this work, we study the Lifshitz transitions in the mag¬ 
netic states of the periodic Anderson model with finite U 
by applying the variational Monte Carlo method.In this 
method, we do not introduce approximations in evaluating 
physical quantities, while we assume variational wavefunc¬ 
tions as in the slave-boson mean-field and Gutzwiller approx¬ 
imation methods. We investigate both the antiferromagnetic 
and ferromagnetic states on an equal footing by varying the 
electron filling. In particular, we analyze the physical quanti¬ 
ties and energy gain at the Lifshitz transitions to determine the 
characteristics of the transitions. Preliminary results on the to¬ 
tal energy for both magnetic cases and the ordered moment in 
the antiferromagnetic case have been reported in Ref. 40. 

This paper is organized as follows. In Sect. 2, we explain 
the periodic Anderson model and the variational wavefunc¬ 
tions in this study. In Sect. 3, we show the calculated results 
for an antiferromagnetic case (Sect. 3.1) and for a ferromag¬ 
netic case (Sect. 3.2). We calculate the energy and physical 
quantities such as the ordered moment and effective mass. We 
also discuss the nature of the phase transitions in the magnetic 
phases with the aid of analyses of the energy components and 
the momentum distribution functions. Then, we discuss the 
Fermi surface structures in the magnetic phases. The last sec¬ 
tion is devoted to a summary. 
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2. Model and Method 

The periodic Anderson model is given by 

^ ^kc\^Ckcr + ^ 

kCT i<T 

+ c^Jka) + U 2^ «/.T«/4^ 

her i 

where and are the creation operators of the conduction 
and / electrons, respectively, with momentum k and spin <j. 
n ficr is the number operator of the / electron with spin cr at site 
i. ek is the kinetic energy of the conduction electron, ef is the 
/-electron level, V is the hybridization matrix element, and U 
is the onsite Coulomb interaction between / electrons. Here, 
we consider only the nearest-neighbor hopping for the con¬ 
duction electrons on a square lattice, and the kinetic energy 
is given by e/c = -2f(cos kx + cos ky), where t is the hopping 
integral and we set the lattice constant as unity. 

We apply the variational Monte Carlo method to the 
model.^^’^® As the variational wavefunction, we consider the 
following Gutzwiller wavefunction: 

I/) - P\^), ( 2 ) 


where 

(3) 


is a projection operator with the variational parameter g. This 
parameter controls the probability of the double occupancy of 
the / electrons on the same site. In the limiting cases, g = 1, 
i.e., P - \ for U - 0 and g - 0, i.e., the double occupancy 
is prohibited for t/ — > oo. For a finite U as in this study, we 
have to determine g between zero and unity to minimize the 
energy. |/) is the one-electron part of the wavefunction. In the 
present study, we choose the one-electron part as the ground 
state of a mean-field-type effective Hamiltonian. 

For the paramagnetic or ferromagnetic state, i.e., for a uni¬ 
form state, we consider the following effective Hamiltonian: 




kcT 


' ek 

-vfi 
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where Va- is the effective hybridization matrix element and 
ifer is the effective /-level. They are variational parameters. 
For the paramagnetic state, they do not depend on spin cr. 

For the antiferromagnetic state, we consider a state with the 
ordering vector Q = (tt, tt). Then, the effective Hamiltonian is 
given by 
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where k-summation runs over the folded Brillouin zone of the 
antiferromagnetic state and cr in front of the parameters stands 
for + (-) for the up-spin (down-spin) states. The parameters 
with a tilde are variational parameters. Ccq and e/g play roles 
similar to mean fields. In addition, we consider Vq, which de¬ 
scribes the staggered component of the effective hybridization 
matrix element in the antiferromagnetic state. 

For the charge-density-wave state with Q = (n, tt), we can 
also consider a similar effective Hamiltonian, but we find that 
the charge-density-wave state does not become the ground 
state within the parameters that we have investigated. In the 
Kondo lattice model, the possibility of the charge-density- 
wave state has been discussed."^^^"* To discuss this possibil¬ 
ity in the periodic Anderson model, we need to investigate 
a much wider parameter space, e.g., by varying U, since 
the charge-density-wave state is considered to be realized in 
an intermediate coupling regime in the Kondo lattice model. 
Thus, we show results only for the paramagnetic, ferromag¬ 
netic, and antiferromagnetic states in the following. 

To construct |/), we fix the number of electrons per site of 
each spin cr, ng-. In the paramagnetic and antiferromagnetic 
states, = np For the ferromagnetic state, the magnetization 
M = n-f - is a parameter characterizing the state. 

For each state, we evaluate the energy by the Monte Carlo 
method, and optimize the variational parameters that min¬ 
imize the energy. Then, we compare the energies of these 
states with the same electron density n - n-i + ni and de¬ 
termine the ground state. Other physical quantities can also 
be calculated by the Monte Carlo method with the optimized 
variational parameters. 

In this study, we set U -%t and V - t, that is, U is the same 
as the bandwidth of the conduction electrons and V is much 
smaller than the bandwidth. The calculations are carried out 
for an LxL lattice with L - 12. The boundary condition is an- 
tiperiodic for the x-direction and periodic for the y-direction. 

3. Results 

5.7 Around half-filling: n — 1.917 

First, we show the results around half-filling (n = 2). We 
set the number of electrons per site n to 276/12^ = 1.917. 

Figure 1 shows the energy E per site of the antiferromag¬ 
netic (AF) and ferromagnetic (FM) states measured from that 
in the paramagnetic (PM) state Lpm as a function of e/. For 
the ferromagnetic states, we show the results for M - 0.083, 
0.583, and 1. The state with M - 0.083 is the half-metallic 
state for this filling, i.e., M - n-^ - n^ - 1 - {n - 1) - 2 - n. 

In a wide parameter region, we find that the antiferromag¬ 
netic state is the ground state. At e//f > -0.1, the half- 
metallic state with M - 0.083 has the lowest energy, while 
the difference in energy is not visible on this scale. The en¬ 
ergy gain of this weak ferromagnetic state is very small, and 
it may become unstable against the paramagnetic state when 
we improve the variational wavefunction. Thus, we simply 
ignore this ferromagnetic state here and concentrate on the 
antiferromagnetic state. In the antiferromagnetic state, there 
is a bend in the energy at ey/f ^ -1.3. The discontinuity in 
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Fig. 1. (Color online) Energy as a function of ef measured from that of the 
pai'amagnetic state FpM for the antiferromagnetic (AF) state (open squares) 
and for the ferromagnetic (FM) states with M = 0.083 (solid squares), M = 
0.583 (circles), and M = 1 (triangles). Ult = 8, Vjt = 1, and n = 1.917. 
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Fig. 2. (Color online) (a) Antiferromagnetic moment Map and (b) elfective 
mass as functions of 6/ for VIt = 8, Vlt = 1, and n = 1.917. The vertical 
lines denote the phase boundaries. We draw the Fermi surface in each phase 
in (a): only the lower hybridized band is occupied in the lightly shaded areas 
and both the hybridized bands are occupied in the darkly shaded areas. 


the first derivative of the energy indicates a first-order phase 
transition. 

In Fig. 2(a), we show the antiferromagnetic moment as a 
function of ej. The antiferromagnetic moment is defined as 

^ “ ”4)’ (6) 

i 

where V = is the number of lattice sites, r, is the posi¬ 
tion of site i, riia- is the number operator of the electrons with 


spin cr at site i, and (■ ■ ■) denotes the expectation value. By 
decreasing ef. Map develops from zero around ej ^0.1. This 
seems to be a continuous phase transition, although we can¬ 
not discriminate it from a weak first-order transition in the 
present numerical calculation. By decreasing e/ further, we 
find ajump in Map at e//f ^ -1.3. This is a first-order phase 
transition as is already recognized from the energy (Fig. 1). 

Here, we call the antiferromagnetic phase with smaller Map 
( ey/f > -1.3) AFl and that with larger Map (ey/f :$ -1.3) 
AF2. For each phase, we can draw the Fermi surface by using 
the obtained variational parameters in the one-electron part 
[see insets in Fig. 2(a)]. We will discuss these Fermi surface 
structures later. 

In Fig. 2(b), we show the effective mass m* defined by the 
jump An(ke) in the momentum distribution function n{k) (see 
Fig. 5) at the Fermi momentum kp: 


m An{kp) ’ 

where m is the bare mass. Here, An{kp) is defined as the jump 
in the momentum distribution function along {n, 0)-(7r, n) for 
the paramagnetic state and along (0,0)-(;7r, 0) for the antiferro¬ 
magnetic states. In the paramagnetic state, m* increases as ey 
decreases, since the number of / electrons increases and cor¬ 
relation effects become stronger. At the PM-AFl phase transi¬ 
tion, m* does not change significantly since it is a continuous 
transition. In the AFl state, m* continues to increase except 
around the PM-AFl phase boundary. On the other hand, in the 
AF2 state, the effective mass becomes lighter since the mag¬ 
netic moment develops sufficiently and the correlation effects 
become weak. 

Note that the symmetry is the same between the AFl and 
AF2 states. To determine what characterizes the AF1-AF2 
transition, we decompose the energy into four terms; the ki¬ 
netic energy of the conduction electrons, 

^ Z (8) 

ko- 

the site energy of the / electrons, 

Eef = eftif, (9) 

the hybridization energy, 

= (10) 

her 

and the Coulomb interaction, 

U 

^ ( 11 ) 

i 

where ny is the expectation value of the number of / electrons 
per site. Figure 3 shows the decomposed terms as functions of 
ey. At the PM-AFl transition, these terms change smoothly. 
At the transition from AFl to AF2, the gain in the hybridiza¬ 
tion Ev decreases, while the gains in E, and E^^ increase. This 
indicates that the conduction and / electrons are relatively de¬ 
coupled in the AF2 state. The change in E^ at the AF1-AF2 
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Eflt Ej-lt 

Fig. 3. (Color online) Components of energy as functions of 6y: (a) kinetic 
energy of the conduction electrons, (b) site energy of the / electrons = 
€fnf measured from ey, (c) energy of the hybridization, and (d) energy of the 
Coulomb interaction. U/t = 8, V/t = 1, and n = 1.917. 



e/li 

Fig. 4. (Color online) Contributions of the conduction and / electrons to 
the electron number and to the antifen'omagnetic moment as functions of ey. 
(a) Numbers of conduction electrons, ric (squares), and / electrons, n/ (cir¬ 
cles), per site, (b) Total antiferromagnetic moment, Maf (open squares), the 
antiferromagnetic moment of the conduction electrons, M^af (solid squares), 
and the antiferromagnetic moment of the / electrons, M/^af (circles). Ult = 
8, V/t = 1, andrt = 1.917. 


transition is small in comparison with the other terms. 

In Fig. 4(a), we show the ef dependences of the occupan¬ 
cies of the conduction and / electrons, is the expectation 
value of the number of conduction electrons per site. The /- 
electron number nj increases as the / level decreases, and in 



(0,0) (7t,0) (7t,Jt) (0,0) (0,0) (7t,0) (7t,7t) (0,0) 


Fig. 5. (Color online) Momentum distribution functions n{k) (open 
squares), tic/k) (solid squares), and nf{k) (circles) for (a) ef/t = 0.5 (PM), (b) 
ef/t = -0.5 (AFl), and (c) ef/t = -2 (AF2). U/t = 8, V/? = 1, and n = 1.917. 
Owing to the antiperiodic boundary condition for the A:-direction, we shift kx 
by tt/L, e.g., {n, tt) in the hgures actually means the point {n — tt/L, n). 


the AF2 state, it almost reaches unity. In Fig. 4(b), we show 
the antiferromagnetic moments of the conduction electrons, 
Mcaf, and of the / electrons, M/af, as functions of e/. They 
are dehned as 

M,af= (12) 

i 

MfAP = <«/4 - «/4 >> ( 13 ) 

i 

where nda- is the number operator of the conduction electron 
at site i with spin cr. In the periodic Anderson model, the con¬ 
duction and / electrons tend to have spins that are opposite 
to each other at the same site. Thus, Mcap and M/af have op¬ 
posite signs. The total magnetic moment Maf is mainly com¬ 
posed of the / component. In the AF2 state, M/af is near to 
unity. This also indicates that the / electrons are almost local¬ 
ized in the AF2 state. 

The momentum distribution functions in each phase are 
shown in Fig. 5. For the paramagnetic or ferromagnetic state, 
the momentum distribution functions are defined as 


ricAk) = {cl^Cka-), (14) 

Hfaik) = {fljk.), (15) 

ria-ik) = rico-ik) + nfa-{k). (16) 

For the antiferromagnetic state; 

f^ccrik) — ^k+Qo-^^ ( 12 ) 

nfcr{k) — {f/j-fkcr fk+Qa-f^+Q°''^’ ( 1 ^) 

tia-ik) = ticAk) + nfa-{k). (19) 


They do not depend on the spin cr in the paramagnetic and 
antiferromagnetic states; nc-[{k) - ric^ik) - ndk), nfiik) - 
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PM: large FS PM: small FS + AF: localized/ 
itinerant/ (fictitious) 


folding of Brillouin zone 



AFl AF2 


Fig. 6. Fermi surface (FS) structure in each phase obtained for Ult = 8, 
Vjt = 1, and n = 1.917. Only the lower hybridized band is occupied in 
the lightly shaded areas and both the hybridized bands are occupied in the 
darkly shaded areas. We also draw the Fermi surface structure for a fictitious 
localized /-electron state with a small Fermi surface composed only of the 
conduction electrons. 


rifiik) = n/(A:), and nd*) = «i(*) = «(*)■ 

In Fig. 5, we recognize most of the Fermi momenta on the 
symmetry axes by the clear jumps in n{k) even in the finite- 
size lattice in the present study. While n{k) should also have 
jumps around (7r/2,7r/2) in the AF2 state (see Fig. 6), we 
could not detect them in the lattice with the present size. In 
the PM and AFl states, the jumps in the total momentum dis¬ 
tribution function n{k) are mainly composed of the / contri¬ 
bution njik). On the other hand, in the AF2 state, the jumps 
are mainly due to the conduction-electron contribution ndk). 
In the AF2 state, nf(k) is almost flat, that is, the / electrons 
are nearly localized in the real space. 

In Fig. 6, we show the Fermi surface structure in each state. 
In the paramagnetic state, there is a small hole pocket around 
{7r,7T) since it is near half-filling. In the present theory, the 
paramagnetic state is always regarded as an itinerant / state, 
that is, the /-electron state contributes to the volume of the 
Fermi surface. In the AFl state, we obtain a hole pocket cen¬ 
tered at (0,0). This Fermi surface can be obtained by simply 
folding the paramagnetic Fermi surface. Thus, the AFl state 
is naturally connected to the paramagnetic state, and in this 
sense, it is regarded as an itinerant / state. In the AF2 state, 
the Fermi surface is different from that in the AFl state. Thus, 
we can discriminate these antiferromagnetic states on the ba¬ 
sis of the Fermi surface structures, while the symmetries of 
these states are the same. 

The Fermi surface in the AF2 state can be obtained by con¬ 
sidering a fictitious small Fermi surface state. If each site has 
one perfectly localized / electron decoupled from the conduc¬ 
tion electrons and these / electrons order antiferromagneti- 
cally, then we obtain a small Fermi surface composed only of 
the conduction electrons with filling n -1. By combining these 
conduction- and /-electron states, we obtain the same Fermi 
surface as in the AF2 state. This means that the AF2 state can 



Fig. 7. (Color online) Energy as functions of 6/ measured from that of the 
paramagnetic state E^m for the antiferromagnetic state (open squares) and 
the ferromagnetic states with M = 0.25 (solid squares), M = 0.5 (circles), 
M = 0.75 (triangles), and M = \ (diamonds). U/t = 8, V/t = 1, and n = 1.5. 


be interpreted as a localized / state. Note, however, that the 
conduction and / electrons are not completely decoupled. 

The AF1-AF2 transition is of first order, since the AF2 
Fermi surface cannot be obtained by continuously deform¬ 
ing that in the AFl state. The PM-AFl transition can become 
of first order in general, while in the present calculation, it is 
continuous. 

In CeRhi_iCoj^In 5 , there are two antiferromagnetic phases 
as in the present theory. The change in the Fermi surface 
between the antiferromagnetic phases is observed by the de 
Haas-van Alphen measurement.'^ The variation in the effec¬ 
tive mass deduced from the de Haas-van Alphen measure¬ 
ment as a function of x is similar to that shown in Fig. 2(b). 
While the transition between the antiferromagnetic phases in 
CeRhi_j^Coj^In 5 is a commensurate-incommensurate transi¬ 
tion, the present theory should have some relevance to this 
material, for example, the mechanism of the change in the ef¬ 
fective mass. 

3.2 Far away from half-filling: n — 1.5 

Next, we show the results for n = 1.5. We expect that the 
antiferromagnetism becomes weak for such a case far away 
from half-filling and there is a chance of stabilizing a ferro¬ 
magnetic state. Figure 7 shows the energy as functions of e/ 
for n - 1.5. In contrast to the case around half-filling, the 
ferromagnetic state has a lower energy than the antiferromag¬ 
netic state in a wide parameter region. We note that while the 
antiferromagnetic state has the lowest energy at e//f = -3 in 
Fig. 7, ferromagnetic states with M ^ 0.9 (not shown) have 
lower energy there. 

To determine the magnetization M for each e/, we calculate 
the energy as a function of M. In Fig. 8, we show results for 
efit - -3 and -1 as examples. For ejlt - -3 [Fig. 8(a)], 
the energy becomes minimum at M ^ 0.89. For e//f = -1 
[Fig. 8(b)], the energy becomes minimum at M = 0.5. The 
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Fig. 8. (Color online) Energy as a function of magnetization M: (a) efjt = 
—3 and (b) efjt = -1. Ujt = 8, Vjt = 1, and n = 1.5. 


state with M - 0.5 is the half-metallic state for this filling. 
We find a cusp in the energy at the minimum point M - 0.5. 
This indicates a gap in the spin excitation for the half-metallic 
state, since the magnetic susceptibility x is given by - 

X^^ and a cusp in E{M) results in;^f = 0. This gap originates 
from the hybridization gap between the up-spin bands. 

In Fig. 9(a), we show the magnetization as a function of 
ef. By decreasing ey, M gradually develops from zero around 
ef ^ 0.4. For -1.7 < e/ < -0.8, we obtain the half-metallic 
state, M = 0.5. The magnetization is flat in this region. By 
decreasing e/ further, M increases again and asymptotically 
reaches unity. In the following, we call the low-magnetization 
state (M < 0.5) FMO, the half-metallic state (M - 0.5) FMl, 
and the high-magnetization state (M > 0.5) FM2. These fer¬ 
romagnetic states have the same symmetry, but we can dis¬ 
criminate them on the basis of the Fermi surface structures as 
shown in Fig. 9(a). We will discuss the details of the Fermi 
surface structures later. 

In Fig. 9(b), we show the e/ dependence of the effective 
mass for each spin state. Here, the effective mass is defined 
along (jT,0)-(n,7T) for the PM, FMO, and FMl phases and 
along (0, 0)-(7T, 0) for the FM2 phase. Note that in the half- 
metallic phase FMl, there is no Fermi surface for the up-spin 
state and we cannot define the effective mass for the up-spin 
electrons. In the PM, FMO, and FMl states, the effective mass 
increases as e/ decreases except at the FMO-FMl boundary. In 
the FM2 state, m* decreases as e/ decreases since the ordered 
moment becomes large. 

In Fig. 10, we show the components of the energy. The 
changes in these components at the phase boundaries are 
weak except for the FM1-FM2 transition. At the transition 
from FMl to FM2, the gain in the hybridization energy Ey is 
reduced, while the gains in the kinetic energy E, of the con¬ 
duction electrons and the site energy of the / electrons 
increase. This indicates that the conduction and / electrons 
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Fig. 9. (Color online) (a) Magnetization M and (b) effective mass for up- 
spin (squares) and down-spin (circles) states as functions of €f for U/t = 8, 
V/t = 1, and n = 1.5. The vertical lines denote the phase boundaries. In the 
FMl state, the Fermi surface is absent for the up-spin state and we cannot 
define the effective mass for it. We draw the Fermi surface in each phase in 
(a): only the lower hybridized band is occupied in the lightly shaded areas 
and both the hybridized bands are occupied in the darkly shaded area. 




Sflt Sjlt 

Fig. 10. (Color online) Components of energy as functions of ef. (a) ki¬ 
netic energy of the conduction electrons, (b) site energy of the / electrons 
= efUf measured from 6y, (c) energy of the hybridization, and (d) energy 
of the Coulomb interaction. Ult = 8, Vjt = 1, and n = 1.5. 


are nearly decoupled in the FM2 phase, as in the AF2 phase 
of « = 1.917. The change in Fj/ at the FM1-FM2 transition is 
smaller than those in the other terms. 

In Fig. 11(a), we show the e/ dependences of and rif. 
nf increases as ef decreases and reaches almost unity in the 
FM2 phase. In Fig. 11(b), we show the magnetization of the 
conduction and / electrons. Me and M/, respectively, and the 
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Fig. 11. (Color online) Contributions of the conduction and / electrons to 
the electron number and to the magnetization as functions of 6y. (a) Numbers 
of conduction electrons, ric (squares), and / electrons, Uf (circles), per site, 
(b) Total magnetization M (open squares), the magnetization of the conduc¬ 
tion electrons. Me (solid squares), and the magnetization of the / electrons, 
Mf (circles). Ult = 8, Vjt = 1, and n = 1.5. 


total magnetization M. Me and Mf are given by 

A/c = ^ ^<ndT - «c4>, (20) 

i 

Mf = ^ “ "/d)- (21) 

i 

At most data points, Me and Mf have opposite signs. Al¬ 
though at some points, Me and Mf have the same sign, the 
absolute values of Me are very small there. The /-electron 
contribution Mf dominates the total magnetization M, and Mf 
is nearly unity in the AF2 phase. 

In actual situations, we should take different values of the 
g-factors for the conduction and / electrons. Thus, the total 
magnetization is not proportional lo M - Me + Mf. However, 
Me is small and the overall features in the total magnetization 
will not change, e.g., the magnetization will remain almost 
flat in the FMl phase. 

Figure 12 shows the momentum distribution functions in 
each phase. In the PM phase, they do not depend on spin. In 
the FMO phase, the number of up-spin electrons increases and 
the hole Fermi surface around in, n) shrinks for the up-spin 
state. For the down-spin state, the hole Fermi surface should 
become larger, but, owing to the small magnetization and the 
small lattice size in the present study, we cannot detect the 
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Fig. 12. (Color online) Momentum distribution functions na-{k) (open 
squares), rico-ik) (solid squares), and nfo-{k) (circles) for (a),(b) efjt = 1 
(PM), (c),(d) €flt = -0.4 (FMO), (e),(f) eflt = -1.5 (FMl), and (g),(h) 
efjt = -3 (FM2). The left (right) panels show those of up-spin (down-spin) 
states. Ult = 8, Vjt = 1, and n = 1.5. Owing to the antiperiodic boundary 
condition for the .jc-direction, we shift kx by nJL, e.g., (n,!:) in the figures 
actually means {n — njL,-!:). 


change. The number of / electrons is larger than that in the 
PM state, and the contribution of the / electrons increases, 
particularly around the Fermi momenta. In the FMl phase, the 
Fermi surface for the up-spin state disappears as is recognized 
from the absence of jumps in nf{k). In the FM2 phase, the 
jumps in no-(k) at the Fermi momenta are mainly composed 
of nea-(k). Hferik) is nearly flat and the / electrons are almost 
localized in the real space. 

In Fig. 13, we show the Fermi surface structure in each 
state. The Fermi surface in the PM state is what is called 
a large Fermi surface with the /-electron contribution. In 
the FMO phase, the hole Fermi surface of the up-spin state 
shrinks, and in the FMl phase, it disappears. In the FM2 
phase, the up-spin electrons partially occupy the upper band, 
and as a result, the Fermi surface structures for up- and down- 
spin states become similar to each other. 

This Fermi surface in the FM2 state can be understood from 
a localized / picture. The Fermi surface in the FM2 state is 
approximately decomposed into a fictitious localized / state 
with complete polarization and a paramagnetic small Fermi 
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Fig. 13. Fermi surface (FS) structure in each phase obtained for U/t = 
8, V/t = 1, and n = 1.5. Only the lower hybridized band is occupied in 
the lightly shaded areas and both the hybridized bands are occupied in the 
darkly shaded area. The Fermi surface for FMO is obtained at efjt = -0.75 
and that for FM2 is obtained at efjt = —3. In the other phases, the Fermi 
surface does not change with €f. We also draw the Fermi surface structure for 
a fictitious localized /-electron state with a small Fermi surface composed of 
the conduction electrons. 


surface of the conduction electrons with filling ric - n - 
Thus, the FM2 state is regarded as a localized / state. 

In the present calculation, the FMO-FMl transition is of 
first order and the other transitions are continuous. However, 
in general, it is possible for each of them to occur through ei¬ 
ther a first-order transition or a continuous transition, since 
each Fermi surface can be continuously deformed into the 
others. 

In UGe 2 under pressure, there are two ferromagnetic 
phases, which probably correspond to FMl and FM2 in this 
study. We are uncertain whether the FMO phase can be elimi¬ 
nated by tuning the parameters. If we ignore the weak magne¬ 
tization in the FMO state, the overall behaviors of the magne¬ 
tization and the effective mass in Fig. 9 are similar to those as 
functions of pressure in UGe 2 .^^’^’*° In addition, the Fermi 
surface reconstructions at the phase transitions are also ob¬ 
served in the de Haas-van Alphen measurements.^’^ ' Thus, 
we expect that the FM1-FM2 transition in UGe 2 is a Lifshitz 
transition corresponding to the present theory. 

4. Summary 

By applying the variational Monte Carlo method, we have 
investigated both the antiferromagnetic and ferromagnetic 
states of the periodic Anderson model with finite U on an 
equal footing. We have found the antiferromagnetic states 
(AFl and AF2) around half-filling (n = 1.917) and the fer¬ 
romagnetic states (FMO, FMl, and FM2) for a case far away 
from half-filling (n = 1.5). 

The weak magnetic states, AFl and FMO, are naturally con¬ 
nected to the paramagnetic state with a large Fermi surface. 
On the other hand, the large ordered-moment states, AF2 and 
FM2, can be regarded as localized / states with a small Fermi 
surface. This has been confirmed from the behavior of several 
quantities: the effective mass, the energy components such 
as the hybridization energy, the momentum distribution func¬ 
tions, and the Fermi surface structure. We have also found a 


half-metallic state FMl between FMO and FM2 for n - 1.5. 

These magnetic phases are characterized by the Fermi sur¬ 
face structure, and the transitions between them are Lifshitz 
transitions without symmetry breaking. This is consistent 
with the previous studies with {/ — > oo that separately dis¬ 
cussed the antiferromagnetic and ferromagnetic cases. While 
we have not found a feature peculiar to a finite-1/ case, it gives 
justification for the use of t/ —> oo in related theories. 

In the present study, by carefully analyzing several quanti¬ 
ties, we have reached a unified picture of the Lifshitz transi¬ 
tions for both the antiferromagnetic and ferromagnetic cases. 
In particular, we have clearly shown that both the transi¬ 
tions to the large ordered-moment states, AF2 and FM2, are 
itinerant-localized transitions of the / electrons. 

However, in the present theory, we could not obtain a large 
effective mass, since the large ordered-moment states appear 
before the effective mass is enhanced substantially. To attain 
a coherent understanding of the heavy-fermion state and its 
magnetic order, we need further breakthroughs, such as im¬ 
proving the wavefunction and/or revising the model. These 
are important future problems. 
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